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ABSTRACT 

We calculate the two-loop 10 x 10 anomalous dimension matrix 0{ag) involving 
current-current operators, QCD penguin operators, and electroweak penguin opera- 
tors especially relevant for AS" = 1 weak non-leptonic decays, but also important for 
AB = 1 decays. The calculation is performed in two schemes for 75: the dimensional 
regularization scheme with anticommuting 75 (NDR), and in the 't Hooft-Veltman 
scheme. We demonstrate how a direct calculation of diagrams involving 75 in closed 
fermion loops can be avoided thus allowing a consistent calculation in the NDR 
scheme. The compatibility of the results obtained in the two schemes considered 
is verified and the properties of the resulting matrices are discussed. The two-loop 
corrections are found to be substantial. The two-loop anomalous dimension ma- 
trix 0{aas), required for a consistent inclusion of electroweak penguin operators, is 
presented in a subsequent publication. 
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1 Introduction 



Effective low energy Hamiltonians for non-leptonic weak decays of hadrons are usually 
written as linear combinations of four-quark operators. The coefficients of these operators, 
the Wilson coefficient functions, can be calculated in the renormalization group improved 
perturbation theory, as long as the normalization scale fi is not too low. The size of 
these coefficients depends on /i through the QCD effective coupling constant and on the 
anomalous dimensions of the operators in question. Since the operators generally mix 
under renormalization, one deals with anomalous dimension matrices rather than with 
single anomalous dimensions. 

In the case of the AS = 1 Hamiltonian, relevant for instance for the AI = 1/2 rule 
and the ratio e'/e, there are ten operators Qi, i = 1, ... ,10 which have to be considered. 
Consequently one deals with a 10 x 10 anomalous dimension matrix. The ten operators in 
question can be divided into three classes: 

• current-current operators Qi, i = 1,2 originating in the usual W-exchange and 
subsequent QCD corrections, 

• QCD penguin operators Qi, i = 3, . . . ,6 originating in QCD penguin diagrams, and 

• electroweak penguin operators Qi, i = 7, ... ,10 originating in the photon penguin 
diagrams. 



Explicit expressions for these operators are given in eq. ( p.l| ). 



Because of the presence of electroweak penguin operators a consistent analysis must 
involve anomalous dimensions resulting from both strong and electromagnetic interactions. 
Working to first order in a but to all orders in the following anomalous dimension matrix 
is needed for the leading and next-to-leading logarithmic approximation for the Wilson 
coefficient functions, 

A, = ^ ^(0) , ^ X,(0) , ^(1) + ^(1) (I I) 

^ Att^' 4:71^' (47r)2 ^(47r)2^^^ ' ^ 

In the leading logarithmic approximation only the one-loop matrices 7^^°^ and 7^°^ have 
to be considered. Inclusion of next-to-leading corrections requires the evaluation of the 
two-loop matrices 7s^^^ and In addition some one-loop finite terms are needed. 

What is known in the literature are the matrices ^^^^ and 7^°-* 10^1 111 ^ind a 2 x 2 



submatrix of j^^^ involving the current-current operators Qi and Q2 [12, 13|. The purpose 



of the present and a subsequent paper is to complete the evaluation of the matrix j^^^ and 
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to calculate 7£-'. Together with certain one-loop finite terms this will allow the extension 
of the phenomenology of non-leptonic AS = 1 transitions beyond leading logarithms. The 
results for the 6x6 submatrix of ■y^^^ involving current-current and QCD penguin operators 
together with some phenomenological implications have already been presented by us in ref. 



1^ . In the present paper we give the details of these two-loop calculations and generalize 
them to the full matrix 7^^^^. The matrix 7^^^^^ is considered in a subsequent publication 



| T5| . In particular, we present the results for an arbitrary number of colours, A^, which is 



useful for the applications of the expansion. The phenomenological implications of 



these new results will be discussed in ref. \16 



The calculation of the matrices 7^^^'' and 7^^^'' is very tedious. It involves a large num- 
ber of two-loop diagrams which make the use of an algebraic computer program almost 
mandatory Moreover one has to deal with several subtleties which are absent in the 
calculations of ■y^^^ and 7^°^ First of all, the two-loop anomalous dimensions depend on 
the renormalization scheme for operators and the treatment of 75 in D 7^ 4 space-time 



dimensions. It is known for instance on the basis of ref. |13|, that the 2x2 submatrix of 
y^^^ involving Qi and Q2 calculated in the NDR scheme (naive dimensional regularization 
with anticommuting 75) differs substantially from the one obtained in the HV scheme (non 
anticommuting 75 |]T^, |2D[). The same feature is found for the full matrices y^^^ and 
The important point is that the difference between 7*^^^ calculated in two different schemes 
is on general grounds entirely given in terms of 7*^°-' and the finite parts of one-loop dia- 
grams calculated in the two schemes in question. This relation is very useful for checking 
the compatibility of two-loop calculations performed in different renormalization schemes 
and plays an important role in demonstrating the scheme independence of physical quan- 



tities. It is given in eq. (|2.15|) . Another subtle point is the dependence of the two-loop 



anomalous dimensions obtained in the NDR scheme on whether a given operator has been 
put in colour singlet or colour non-singlet form such as 

Q2 = {SaUa)Y_^{Updf3)^_^ OI Q2 = {Sadp)y_j^{Ui3Ua)y_j^ ■ (1-2) 

In D = 4 dimensions Q2 and Q2 would be equivalent by means of Fierz transformation. 
Their one-loop anomalous dimensions are the same. 

Again the two-loop anomalous dimensions for these two forms of operators can on 
general grounds be related to each other by calculating one-loop diagrams. This feature 
turns out to be fortunate in that it offers one way to deal with dangerous closed fermion 



^^The algebraic computer program for the 7-algebra TRACER, written in MATHEMATICA IT7|, which 
we extensively used during the calculation, has been developed by two of us [ |T8[ and is available from the 
authors. 
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loops involving odd numbers of 75 in two-loop diagrams. The latter are known to be 
ambiguous in the case of the NDR scheme. Indeed, as we will demonstrate explicitly in 
section 5 and in ref. [IT3I, it is possible by using different forms of the operators and the 



relation mentioned above, to express all diagrams with closed fermion loops in terms of 
diagrams not containing such loops. In this way, consistent calculations of ■y^^^ and ■y^l^ can 
be performed in the NDR scheme and compared with the one obtained in the HV scheme, 
free of 75 problems. 

The following sections constitute the calculation of 7^^^^ and a detailed discussion of the 



subtle points mentioned above. In ref. |15|, the methods developed here are generalized to 
include photon exchanges and are used to calculate 7^''. 

Our paper is organized as follows: In section 2, we give a list of the ten operators in 
question, and we classify the one and two-loop diagrams into current-current and penguin 
diagrams. We also discuss the basic formalism necessary for the calculation of In 
section 3, we recall the matrix and we calculate the finite terms of one-loop diagrams 
in the NDR and HV scheme. In section 4, the calculations and results for two-loop current- 
current diagrams are presented. In section 5, an analogous presentation is given for two- 
loop penguin diagrams. In section 6, we combine the results of the previous sections 
to obtain 7g('^) in NDR and HV schemes. We discuss various properties of this matrix, 
in particular its Large- iV limit. Section 7 contains a brief summary of our paper. In 
appendices A and B explicit expressions for the elements of the 10 x 10 matrices 7g^°^ and 
T-g^^-* for an arbitrary number of colours (A^) and flavours (/) are given. In appendix C the 
corresponding results for 7^^^ in the case of = 3 are presented. 

2 General Formalism 
2.1 Operators 

The ten operators considered in this paper are given as follows 

Ql = (SaM/3)v_A (^/3C?a)v-A ' 

Q2 = (sm)v_a Mv-A ' 

Q3 = (sC?)v-A i^lh-A ' 
Q 

Q5 = (S(i)v-A J2 (^?)v+A ' 
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Qe = 


is, 




q 


Qi = 


3 




-A E e<? (^?)v+A ' 
<? 


Qs = 


3 
2 


[Sadp] 


•v-aX!^9 (^/3?q)v+A 


Q9 = 


3 
2 


{sd)y^ 


~A E ^9 (^g)v-A ' 
9 


Qio = 


3 
2 


[Sadp] 


'v-A E ^9 (^/3?")v-A 
9 



where a, P denote colour indices = 1, . . . , A^) and are quark charges. We have 

omitted colour indices for the colour singlet operators. (V ± A) refers to 7^(1 ± 75). This 
basis closes under QCD and QED renormalization and is complete if external momenta 
and masses are neglected. However, at intermediate stages of the calculation, we have to 
retain operators which vanish on-shell. This will be discussed in detail in section 5.2. 

At one place, it will become useful to study a second basis in which the first two 
operators are replaced by their Fierz conjugates, 

Qi = (^C^)v-A (^«)v-A ' 

Q2 = {Sadi3)y_j^ (M/3Ma)v-A ' (2-2) 

with the remaining operators unchanged. In fact, the latter basis is the one used by Oilman 
and Wise 0. We prefer however to put Q2 in the colour singlet form as in eq. ( |2.1|) , because 
it is this form in which the operator Q2 enters the tree level Hamiltonian. Let us finally 
recall that the Fierz conjugates of the (V — A) ® (V + A) operators Qq and Qg are given 
by 

Qg = (^l^r) (^R^l) , 

9 

Qs = -12^e,(sLgR)(gR4) , (2.3) 

9 

with similar expressions for Q5 and Qj. Here qn^i = ^{l±'y^)q. The Fierz conjugates of the 
(y —A)®(y —A) penguin operators Q3, Q4, Qg and Qio, to be denoted by Qi, i = 3,4, 9, 10, 
can be found in analogy to (|2.2|). 

2.2 Classification of Diagrams 

In order to calculate the anomalous dimension matrices ■y^^^ and one has to insert 
the operators of eq. ( f^.lD in appropriate four-point functions and extract 1/e divergences. 
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The precise relation between the 1/e divergences in one- and two-loop diagrams and the 
one- and two-loop anomalous dimension matrices will be given in the following subsection. 
Here, let us only recall that insertion of any of the operators of eq. ( ^.1| ) into the diagrams 
discussed below results in a linear combination of the operators Qi. A row in the anomalous 
dimension matrix corresponding to the inserted operator can then be obtained from the 
coefficients in the linear combination in question. 

There are three basic ways a given operator can be inserted into a four-point function. 
They are shown in fig. 1, where the cross denotes the interaction described by a current 
in a given operator of eq. and the wavy line denotes a gluon. We shall refer to the 

insertions of fig. 0(a) as "current-current" insertions. The insertions of fig. |l|(b) and (c) 
will then be called "penguin insertions" of type 1 and type 2 respectively. 

The complete list of the diagrams necessary for one- and two-loop calculations is given 
in figs. 2-7. At the one-loop level one has three current-current diagrams, fig. 2, to be 
denoted as in ref. |jl3| by Di - and one penguin diagram of each type, fig. 3, to be 
denoted by Pq^^ and Pq'^ for type 1 and type 2 insertion respectively. At the two-loop level 
there are 28 current-current diagrams shown in fig. 4, to be denoted by D^^ - D^i, and 14 
penguin diagrams of each type shown in fig. 5, to be denoted by Pi^^ - Pil ^ and pI^^ - Pif 
respectively. The diagrams obtained from the ones given here by left-right refiections have 
not been shown. In the case of current-current diagrams also up-down refiections have 
to be considered. In fig. 6, we show two penguin diagrams which have no 1/e divergence 
and hence do not contribute to the anomalous dimensions, and in fig. 7, some examples of 
penguin diagrams which vanish identically in dimensional regularization are given. 



2.3 Basic Formulae for the Anomalous Dimensions 

The anomalous dimensions of the operators Qi, calculated in the MS scheme, are obtained 
from the 1/e divergences of the diagrams of figs. 2-5 and from the 1/e divergence of the 
quark wave function renormalization. Let us denote by Q a column vector composed of 
the operators Qi. Then 

7((7) = Z-^fi^Z, with = ZQ, (2.4) 

Ct/i 

where Q^ stands for bare operators. Working in D = A — 2e dimensions, we can expand 
Z in inverse powers of e as follows 

OO -j^ 

Z = ^ + E7^Z,ig). (2.5) 

k=i ^ 
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where g is the QCD couphng constant. Inserting ( |2.5| ) into ( |2.4| ) one derives a useful 
formula 



l{9) = -'29' 



(2.6) 



Let us next denote by r^'^\Q) and r^^\Q^) the renormalized and the bare four-quark 
Green functions with operator Q insertions. Strictly speaking T^'^^Q) and T^^\Q^) are 
matrices, because the insertion of a single operator in a given diagram results in a linear 
combination of operators. 

At the one-loop level Fg'' is obtained by evaluating the diagrams of figs. 2 and 3. At 
the two-loop level it is found by evaluating the diagrams of figs. 4 and 5 and subtracting 
the corresponding two-loop counter terms. For the results of tabs. 1 - 4, we have made 
these subtractions diagram by diagram. Next 



F(^)(Q) 



(2.7) 



where is a matrix which represents the renormalization of the four quark fields. In the 
case of pure QCD it is a diagonal matrix with all the diagonal elements being equal to 
Z^*^) , where Z^^^^ is the quark wave function renormalization, defined by 



B 



We expand F^"* and Z^"^^ in inverse powers of e as follows 



where 



and 



00 



i + Y.-kZP{g,Q 



finite . 



k=l 



167r2 



e 



+ 02 



167r 



2^2 



zPig,Q^ 
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E 



(2.8) 

(2.9) 
(2.10) 

(2.11) 
(2.12) 



Demanding F*^^) to be finite, we find Zi, and, using (|2.6| ), the basic formulae for the 
one- and two-loop anomalous dimension matrices. 



(7, 



-2 
-4 



2ai6ij + {bi 
2a2Sij + (62' 



(2.13) 
(2.14) 
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2.4 Renormalization Scheme Dependence of 7g(^) 

The two-loop anomalous dimension matrices depend on the renormalization scheme and 



in particular on the treatment of 75 in D 7^ 4 dimensions. In ref. we have derived a 
relation between (7s^-')a and (78^-*);, calculated in two different renormalization schemes a 
and 6, 



(7: 



(7: 



+ 



Ars,7 



(0) 



2/?oArs 



with Afc 



[rsb 



where /3o is the leading coefficient in the expansion of the /3-function 



Pig) = 

and the {fs)i are defined by 

{Q)^ 



-Po 



167r2 



g-' 



:i67r 



2^2 



(2.15) 



(2.16) 



(2.17) 



Here, Q'-^-' is a tree-level matrix element and {Q)i denotes renormalized one-loop matrix 
elements calculated in the scheme i. The matrices (fs)j are obtained by calculating the 
finite terms in the one-loop diagrams of figs. 2 and 3. 

Relation (|2.15|) is very useful as it allows to test compatibility of the two-loop anomalous 
dimension matrices calculated in two different renormalization schemes. It also plays a 
central role in demonstrating the scheme independence of physical quantities [|1^ . 



2.5 Collection of Useful Results 

Let us recall the values for /5o, cti, and 02, 
/3o 



11 2 
—N--f 
3 3'' 



A = |iV^-yiV/-2C^/^ 



ai 



02 



4 4 2'^ 



(2.18) 
(2.19) 



where Cp = [N"^ — 1)/2A^. Here, is the number of colours and / the number of active 
flavours. All four quantities in eqs. ( |2.18| ) and ( |2.19| ) are common to the NDR and HV 
scheme. Moreover Pq and /?i are gauge independent. The coefficients ai and 02 given here 
correspond to the Feynman gauge. 

In the course of the discussion of the anomalous dimensions of the penguin operators 
Qq and Qs it will be useful to have at hand the anomalous dimension of the mass operator 



^ = ^(0)^,^(1) 
Im fm ^Q^2 ^ I'm 



(IGvr 



2^2 



+ 



(2.20) 
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To two loops, the anomalous dimensions 7^-' and 7^^ are given by 



/ m 



6Cf 



and 



7. 



(1) 



(2.21) 



for both schemes 121 



Finally let us recall that for the anomalous dimension of the weak current, 

+ ■ ■ ■ . 



(0) 9 , (1) 9 
7. = 7. Y^ + 7. 



(167r 



2^2 



f2.22) 



the one-loop coefficient 7j'^ vanishes in both schemes. However, at the two-loop level in 



the HV scheme, 7j is found to be non-vanishing |T^, 



if 




NDR 
HV 



(2.23) 



In the HV scheme the axial current attains a finite renormalization if one subtracts mini- 
mally, as we did, and one gets a mixing between chiral components. This is not so aesthetic 
and on hindsight it would have been more elegant to redefine the currents non-minimally to 
retain standard normalization. We do not proceed along this line but in our computation 
the mixing reappears in the non- vanishing of the two-loop coefficient '^^P ■ However, since 
the one-loop matrix elements f also change correspondingly, the final physical quantities 
are normalized properly and scheme independent in the end. 



3 One— Loop Results 

3.1 Current— Current Contributions to 

The contributions of diagrams Di - D3 of fig. 2 to the matrix 'y^^^ have a very simple 
structure. The mixing between different operators can be divided into five blocks, each 
containing two operators, 

iQk,Qk+i) , A; = 1,3,5,7,9, (3.1) 

with no mixing between different blocks. Moreover, the mixing among {V — A) ® {V — A) 
operators [k = 1, 3, 9) is described by a universal 2x2 matrix. Using eq. ( |2.13| ), one finds 

ii'\QuQ2) 



Q/N 6 
6 -6/iV 



(3.2) 
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7f (Q5,Q6 



with identical results for {Q3,Q4) and {Q9,Qio)- 
Similarly, 

' 6/N 


with identical result for {QjiQ^)- 

Comparing with eq. ( |2.21|) , we note that 

,(0)1 _ L,(o)l 



-6 

6(iV2 - l)/N 



(3.3) 



2 7. 



(0) 



(3.4) 



showing that the one-loop anomalous dimensions of the operators and Qs are twice 
the one-loop anomalous dimension of the mass operator, provided only current-current 
diagrams Di - are taken into account. It should be pointed out that this result is valid 
for an arbitrary number of colours A^, but remains only true in the Large-A^ limit when 



penguin insertions are taken into account |22, 23 



3.2 Penguin Contributions to 

The contributions of the diagrams Pq^"* and Pq^'' of fig- 3 to the matrix 7^^°^ have also a very 
simple structure: 

• The insertions of Qi, Q^, and Q-j into the diagrams Pq^^ and Pq^^ vanish. This 
follows either from colour conservation, or flavour conservation, or finally from the 
Dirac structure, as one can easily verify by inspecting the diagrams of fig. 3. Thus, 



the corresponding rows in 



7 



(0) 



vanish, 



0. 



(3.5) 



The insertion of any of the remaining operators into the diagrams of fig. 3 results 
always into a unique linear combination of the QCD penguin operators Q3 - Qq, 
multiplied by an overall factor characteristic for the inserted operator. Denoting by 

1 



P 



0, 0, 



1 



N' 



1, 0, 0, 0, 



(3.6) 



the row vector in the space (Qi - Qio), the non- vanishing elements of 7^°^ coming 
from the diagrams of fig. 3, are as follows. 



7fHQ4) 

7f (Qio) 



|(n- 
liu- 



P, 
P. 



7f(Q3) 
7f(Q6) 

7f (Q9) 



fP 



_ Hp 

3 ' 



(3.7) 
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where u and d denote the number of up- and down-quark flavours respectively {u + 
d = f). In obtaining ( p.7|) , eq. ( p.l3|) has been used with ai omitted, since this 



contribution has already been included in (|3.2|) and 



Combining the results of eqs. ( |3.2| ), ( |3.3| ), ( |3.5| ), and (p.7|), we find the known complete 
one-loop anomalous dimension matrix For completeness this matrix has been given 
explicitly in appendix A. 

It should be stressed, that these results do not depend on the renormalization scheme 
for Qi and are also valid for Qi. This is no longer the case for the finite terms in the 
diagrams of figs. 2 and 3, to which we now turn our attention. 

3.3 Current— Current Contributions to Afc 



The contribution to the matrix fg defined in eq. ( p.l7| ), resulting from the diagrams Di - 
D3 of fig. 2. 

• depends on the renormalization scheme for Qi, but 

• does not depend on whether Qi, or their Fierz conjugates Qi, are inserted in these 
diagrams. 

The general structure of [fs\cc is the same as in the case of ■y^'^^ presented in section 3.1. 

Since f is generally gauge dependent and dependent on the infrared structure of the 
theory we give here only results for 

Afs = fs^HV - ^'s.NDR (3.8) 

which are free from these dependences. The results for fg in the Landau-gauge have been 



given in ref. [|T4 
We find 



[Af.W„Q.)L= I '^-//^ I, (3.9) 



with identical results for {Q3,Q4) and {Q9,Qio). Similarly, 



|Ar.W5.ge)L= I '^-/^^ ^^'^^^ I. (3,10) 



with identical result for {Qj, Q 



8) 
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3.4 Penguin Contributions to Afg 

Turning our attention to the finite parts of tfie penguin diagrams of fig. 3, it sliould be 
first noted tliat at the one loop level one does not face the problem of the evaluation of 
Tr[757^7jy7A7w] in D ^ 4 dimensions. Consequently, the diagram P^^^ can be calculated in 
the NDR renormalization scheme without any difficulties. Nevertheless, different results 
are obtained for the NDR and HV scheme. 

The reason that no problem with closed fermion loops arises at this level is as follows. 
Although at the intermediate state Tr(7^757j,7A7^) appears it is contracted with g^i, and 
q^^qi, where q is the momentum of the gluon. Thus actually only the trace Tt^'Jt'J^'Jx) has 
to be calculated, which however is zero. 

The structure of the results is similar to the ones given in section 3.2. In addition to the 
above mentioned scheme dependence, in the case of the NDR, but not the HV scheme, the 
results for the insertions of the operators Q2, Q3, Q^, Q9, and Qio depend on whether these 
operators are taken in the basic form of eq. (|2.1| ) or in the Fierz conjugate form Qi. In the 
case of pure QCD corrections the results for Qi, and for the (V — A) ^ (V + A) operators 
Q5 ^ Qs, do not depend on the form used. These properties will play an important role in 
section 5. 

Beginning with the basis of eq. (|2.1| ), we find that the insertions of Qi, Q5, and Q7 
into the diagrams Pq^^ and Pq^^ vanish in both renormalization schemes considered. Thus 
the corresponding rows in [rs,NDR]p and [fs^Hv]p vanish, 

[rs,NDR(Ql)]p = [fs,NDR{Q5)]p = [^s,NDR(<57)]p = 0, (3.11) 

with an identical result for HV scheme. 

The remaining insertions are non-zero but [Afgjp turns out to vanish for Q4, Qq, Qs 
and Qio: 

[Af3(Q4)]p = [Af3(Q6)]p = [Af3(Q8)]p = [Af3(Qio)]p = 0. (3.12) 
For the remaining operators we find then 

[A^s(g2)]p = -ip, 

[Af,(Q3)]p = -IP, (3.13) 
[Af3(Q9)]p = |P, 

where the vector P has been defined in (|3.6| ). 

Since, as stated above, Q2 and Q2 insertions give different results in the NDR scheme. 
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we find that the Ts.ndr calculated in the bases and 



Afs,NDR = fs,NDR[2.2] - fs.NDRplI 



This result will play an important role in section 5. 



2|) differ. One obtains 
/ \ 







V ; 



(3.14) 



4 Current— Current Contributions to the Two— Loop 
Anomalous Dimension Matrix 

4.1 {V - A)0{V - A) Operators 

The generalization of the one-loop matrix of eq. ( p.2| ) to two-loop level can be done using 
the results of ref. |1^, where the eigenvalues of the 2x2 submatrix {Qi,Q2) have been 
calculated in the NDR and HV schemes. Since also at two-loop level the corresponding 
submatrix is symmetric with equal entries on the diagonal, one has 

1 





cc 






cc 





1 



+ 2 7? 



(4.1) 
(4.2) 



with 7±*^^^ given by formulae ( |5.2| ) and (|5.3|) of ref. [|T^ for NDR and HV schemes respec- 
tively, and 7j^^ in (|2.23| ) of the present paper. 

The details of the calculation of the diagrams D4 - D31 of fig. 4, with {V — A)^(y — A) 
insertions, can be found in ref. fl^, where a table of 1/e^ and 1/e singularities in the 
diagrams 1^4 - D31 has been given. In performing these calculations one has to take care of 
evanescent operators which, although vanishing in D = 4, affect the two-loop anomalous 



dimension matrix of physical operators Qi of eq. (|2.1|) . For this reason, the treatment of 
counter diagrams has to be done with care. The discussion of evanescent operators can 
be found in ref. |TB[, and a nice presentation is given also in Here it suffices to give 
only the projections on the space of physical operators. Denoting the evanescent operators 
generally by Eais^-ys, with a, P, 7, 6 being Dirac indices, the relevant projection for diagrams 
involving {V — A) ^ {V — A) operators used in ref. [|I^ is given by 



EafS,^5 (7r(l + 75))^^(7"(l + 75))5„ 



0, 



(4.3) 
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where in the HV scheme 7,- has to be taken in D = A. 

Using eqs. (^4.1|) and (^^) , together with the results of ref. we find in the case of 
the NDR scheme 

^NDR -BndR 
-Bndr ^NDR 



where 



A 



NDR 



22 2 / 57 1 
~y ~ 3iV ^ YiV^ ' 
with identical results for (QsyQ^) and (Qg^Qio)- 
For the HV scheme we find 



19,, 2, 39 



^HV -ShV 

-Bhv ^hv 



where 



44 



HV 



110 14 / 57 1 



3 ' 3' ' 3 3 2 
with identical results for (QsyQ^) and (Qg^Qio)- 



5 



HV 



2 •' at 



(4.4) 
(4.5) 



(4.6) 



(4.7) 



4.2 (y - ^) (g) (y + A) Operators 

The calculation of the {V — A) ^ {V + A) insertions into the diagrams - D31 proceeds 
in analogy with the one for {V — A) ® {V ~ A) operators. However, the projection on the 
space of physical operators takes this time the following form, 



Eaf3,^S (1 - 75)/37 (1 + 75)^0 



0. 



(4.^ 



In order to check Fierz symmetry properties, it is instructive to calculate simultaneously 
the insertions of the (1 + 75) (1 — 75) operators, e.g. Qq of eq. (p.3|). In this case, the 
projection on the space of physical operators takes the form 



EaP,^S (7r(l + 75))^^ (7"(l-75))5„ 



0. 



(4.9) 



The singular terms in the diagrams - D^i with and Qq insertions in the NDR 
and HV scheme are collected in table 0. The two-loop counter terms have been included 
diagram by diagram. Also corrections to counter diagrams resulting from the mixing with 
evanescent operators have been taken into account. Colour factors have been omitted. 



They can be found in ref. ||T3[. The overall normalization is such that after the inclusion 
of colour factors separate contributions to (&2)jj of eq. ( p^.l4D are obtained. In the case 
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of diagrams 29 - 31, we have included in table ^ the colour factors present in the gluon 
self-energy 

F, = ^N--f, F2 = -N--f. (4.10) 

1 3 , ^ 9 9 ^ ^ 

We make the following observations: 

• the 1/e^ singularities are the same in both schemes as expected, 

• the 1/e singularities differ, which implies different two-loop anomalous dimension 
matrices, 

• the Fierz symmetry in diagrams 1^4 - D^i is preserved in both schemes. The Fierz 
conjugate pairs of diagrams are (4,6), (7,9), (10,12), (13,15), (16,21), (17,20), (18,19), 
(22,24), (25,27), and (29,31). 



The remaining diagrams are self-conjugate. In ref. ||13[, where the operators Qi and Q2 
have been considered the Fierz conjugate diagrams with Qi or Q2 insertions had the same 
singular structure. Now, the Fierz symmetry relates the Q5 insertion of a given diagram 
to the Qq insertion of the Fierz-conjugate diagram. As seen in table |l], this relation is 
satisfied for all Fierz conjugate diagrams: Q5 and insertions give the same result. We 
should, however, emphasize that as in ref. ||13|, also here this structure is only obtained 
after the evanescent operators have been properly taken into account. Otherwise, the Fierz 
symmetry in pairs (16,21), (17,20), (18,19) is broken and the results for diagrams 5 and 
23 are modified. Clearly, the fact that the results in table |l| satisfy all these relations, is a 
good check of our calculations. 

Including colour factors, summing all diagrams D4 - D31, and using the formula ( |2.14| ), 
gives the current-current contribution to the two-loop anomalous dimension matrix of 
{V — A) ^ {V + A) operators. For the NDR scheme we find 



7i,NDR(Q5, Qe 



Cndr -Dndr 

-^NDR -^NDR 



(4.11) 



where 



_ 137 22 / 15 1 
Cndr - ' Yn ^ YPP 
71 18 

^NDR = -yiV + 4/--, 



/^NDR = -^^ + f/ + |' (4-12) 

^NDR = 27« + f -4^ + ^, (4.13) 



with identical results for (QjjQs)- Here, 7^-* is the two-loop anomalous dimension of the 
mass operator given in eq. ( |2.21j ). 
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For the HV scheme we find 




(4.14) 



where 
Chv 



44, . 8,,^ 71 10 / 15 1 

— N Nf \ — H 

3 3^ 6 3 N 2 



/^HV^-fiV-^Z + A, (4.15) 



6 3'^ 



HV 



2 7. 



88,. 16,,^ 229 20 / 9 

H Nf \ — H 

"^ 3 3-^ Q 3 N m 



(1) 



with identical results for [Q-j, 



4.3 Properties and Test of Compatibility 

The matrices given in eqs. ( [4.4|) , ([4.6|) , (|4.11| ), and (|4.14|) give the part of the two-loop 
anomalous dimension matrix coming from the current-current diagrams of fig. 4. Let us 
denote this part by 7^ ndr ^"^^ 7i,HV NDR and HV schemes respectively. We 

note, that 



these two matrices differ considerably from each other. In particular all diagonal 
entries in the HV scheme increase as A^^, whereas in the NDR scheme this is only 
the case for the elements (6,6) and (8,8). 

Since 7gg^ and 78g'' grow like A^ for large A^ and 7gg'' and 78g'' grow like A^^ in both 
schemes considered, we find remembering as ~ 1/^, that these two elements ap- 
proach non- vanishing constants in the Large- A^ limit. In the NDR scheme all other 
elements vanish in this limit and this fact will not be changed by including penguin 
diagrams. In the HV scheme all diagonal elements survive the Large- A^ limit. 

JO) ^,(0) 



the simple relation between 7gg\ 7gg^ and 7^-', given in (|3.4| ), is violated at two-loop 



level in both schemes. However, in the NDR scheme it is recovered for large A^ |2^ 



the vanishing elements (5,6) and (7,8) in eq. ( p.3|) become non-zero at the two-loop 
level. 

using the results of this and the previous section, it is interesting to note that the 
relation ( p.l5| ) is satisfied even if only current- current contributions to and 
7s^^-' are taken into account: 
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+ 2PoA[hL ^ (4-16) 
where 

A [^Jcc = [^s,Hv]cc - [^s,NDr]cc (4-17) 

This is easy to understand. Indeed, the two-loop current-current diagrams of fig. 4 
have subdiagrams of current-current type only, and consequently the compatibility of the 
two calculations should be verified within current-current diagrams alone. The fact that 
this is indeed the case is a good check of our results. 



7s,HV 



7s,NDR 



A 



sJcc ' 



7, 



(0) 



5 Penguin Diagram Contributions to the Two-Loop 
Anomalous Dimension Matrix 



5.1 General Structure 

The calculation of the penguin diagram contributions to ■y^^^ can be considerably simplified 
by first analyzing the general structure of these contributions. The two-loop penguin 
diagrams are shown in fig. 5, where two types of insertions of a given operator, b) and 



(1) 



can 



c) of fig. 1, have to be considered. Yet, as we shall show now, the full matrix 7, 
be obtained by calculating just the insertions of the operators Qi, Q2, Q5, and Qe, which 
moreover receive only contributions from the diagrams Pi - P14 , i.e. penguin diagrams 
without closed fermion loops. 

We first note that the anomalous dimensions of the {V — A)®{V — A) penguin operators 
Qi, Q45 Q9, and QiQ can be expressed in terms of the anomalous dimensions of Qi, Q2, 
Qi, and Q2 as follows, 

(5.1) 
(5.2) 

(5.3) 

(5.4) 





p 


= / 




+ 2 

p 






p 


= / 




+ 2 

p 





d 



where the last terms in (EA 



are only due to internal down-quarks. 
Similarly, the anomalous dimensions of the (V — A) ^ (V + A) operators Qj and Qg 



can be expressed in terms of those of Q5 and Qq respectively. 



{u - d/2) 
1 



(5.5) 
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{u-d/2) 



(5.6) 



We next note that Qi, Q2, Q5, and Qq insertions are non-zero only for Pi^^ - 



,(1) 

14 5 



whereas Qi, Q2, Q5, and Qe insertions are non-zero only for P^ 



(2) 



5(2) 
14 • 



In the HV scheme, there is no difficulty in evaluating the diagrams ^ - Pi^P which 



involve closed fermion loops. Consequently, in this scheme, eqs. (|5.1|) - (|5.6|) can directly 



be used to find the matrix 

>(2) 



/s,HV 



. However, in the NDR scheme it is much easier to deal 



with the diagrams Pi'' -PS\ We will now show that relation (|2.15| ) can be used to reduce 
the full calculation to the latter diagrams. 

Let us first observe that eq. (|2.15|) can also be used to relate the two-loop anomalous 
dimensions obtained in the bases (|2.1| ) and ( p.2|) . Since these two bases differ only in the 
two first operators, and the current-current anomalous dimensions of {Qi, Q2) and (Qi, Q2) 
are the same, we find 



7«[(2.2)] 



+ 



(0) 



where 



Af, 



s,NDR 



Afs,7, 



NDR 
HV 



2/?oArs 



with AfsNDR given in eq. ( |3.14D . Eqs. ( p.7[ ) and ( p.8[ ) allow to find \Q 



(5.7) 

(5.8) 
and 



once 



and 



tIndrIQi) 

7iNDR(Q2) 



71ndr(<5i) 
71ndr(<92) 



+ 2P 



are known. In the NDR scheme we then have 

(5.9) 
(5.10) 



9 



+ - 0, 0, 2 - — 



Ar2' 



iiiv+ii, 11+4, 



2iV-^, 0, 0, 0,0 



whereas in the HV scheme, the insertions of Qi and Q2 are equal to the insertions of Qi 
and Q2, respectively. 

Next, we recall the important result of section 3: the finite pieces in one-loop diagrams 
(the matrix did not depend on whether Q5 and Qq or and Qq had been used. Since 
(QsjQe) and {Q5,Qq) form closed sets under the renormalization due to current-current 
diagrams, and the insertion of any of these four operators in one-loop penguin diagrams 
give always linear combinations of - Q& it is evident that 

(5.11) 
(5.12) 





p 






p 
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in both HV and NDR schemes. 

To summarize, we have shown that the full matrix 



can be obtained by calculating 
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just the insertions of the operators Qi, Q21 Qbi and Qq into the diagrams P| - P] 
5.2 Details of the Two— Loop Calculation 

So far, we have restricted our discussion to "on-shell" results. It is well known however, 
that to renormalize operators one has to include mixing with all operators having the same 



quantum numbers. In the literature p6|, |23, |2^, ^ , it has become standard to divide these 



into 3 classes I, IP, and II'', according to 

• class I: gauge invariant operators that do not vanish by virtue of the equations of 
motion, 

• class IP: gauge invariant operators that vanish by virtue of the equations of motion, 
and 

• class II*: gauge non-invariant operators. 

In our case, the situation is even slightly more complicated, since one has to consider also 
mixing with evanescent operators, which vanish when restricted to 4 dimensions. 

The renormalized operators Qi of eq. (|2.1| ), corresponding to the unrenormalized oper- 
ators Qf , which belong to class I, are given by 

6 K 

Qi = H Zr^' Qf + E P^fc Gk + E, + Ni. (5.13) 

j=i k=i 

Here, Ei are the evanescent operators, Ni are in class II'', and Gi are gauge invariant 
operators of dimension 6 that just involve two quark fields. There are K = 3 independent 
operators of this sort (up to total derivatives), having (sd) quantum numbers which arise 
in our computation up to two-loops, 

Gi = is[D^, [D>^,D,]]Yil~l,)d) , 

G2 = {s{D\D,}Yil-l,)d) , (5.14) 
G3 = (rsD^D,D,S^'''\l-^,)d 



where Dn is the covariant derivative and S'^'^^ = 7^7'^7'^ — 'y^'j^'^t^ 



^The operator {s[D'^, D^]7'^(1 — 75) d) doesn't enter at two-loops. 
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The operator Gi is proportional to the penguin operator Qp which occurs at one-loop, 

Qp = -(s7^(l-75)rfD,Fn, (5.15) 
9 

and 

G2, G1 + G3, and Gi-g4-g6 + ^(Q3 + Q5), (5.16) 

all belong to class IP. 

In arbitrary gauges it is known that class I operators generally mix with class II op- 
erators, in particular with those of class 11^. One can simplify the work considerably by 
using the background field gauge |]5D[, for which no such mixing with gauge non-invariant 



operators occurs. In our calculation, we therefore used such a gauge and thus have 
A^j = in eq. (|5.13|) . However this choice only affects the diagrams with triple gluon 
vertices, Pq and Pj of fig. 5. 

In order to include mixing with the operators of class IP, (|5.14| ), we have to calculate 
the operator insertions into four-quark Green functions with arbitrary external momenta. 
This results in additional momentum dependent operators which reflect the presence of 
class IP operators. The various tensor products, leaving out the external quark fields, 
which appear in the calculation are as follows, 

Vi = 7,^(1-75) ® 7^ 
V2 = ^(1-75)® 



Vs = Vi^-, (5.17) 

= ( ^1(1-75) ® ^2+ ^2(1-75) ® ^1) 4 



q2 



= S,xuPlp'2i^-l,) ® 7'4 

gZ 



where pi, p2 are incoming (outgoing) quark momenta and q = pi — P2- The operators Qi 
of eq. ( p^.lD all correspond to the structure Vi. 



The procedure to calculate the diagrams of fig. 5 followed that described in ref. [13|. 
That is, we first classified all independent integrals occuring in the diagrams and their 
singular parts, computed for arbitrary external momenta. Then the contractions with the 
corresponding 7-matrix structures, using the NDR rules and the HV rules, were performed 
separately. For the NDR scheme this was done by "hand" and by using the algebraic 



•^^We did not compute in an arbitrary covariant background field gauge, but restricted ourselves to the 
Feynman background field gauge. 
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computer program TRACER |T^ . For the HV scheme we used TRACER nearly exclusively 
for the penguin diagrams. 

The operator structures V2 and V4 already vanish going on-shell, i.e. using the equations 
of motion, while the on-shell projection of V3 and V5 is given by 

Vs ^ - ^V, and V, V,. (5.18) 

As a check, the on-shell result was also calculated in the ordinary Feynman gauge. Since 
only diagrams Pq and P7 are affected, after using the equations of motion their sum should 
agree with the result in the background field gauge. Up to evanescent contributions this is 
indeed the case. 

A second technicality of the calculation that we like to discuss in somewhat more detail 
is the treatment of 75. 

There are basically two diagrams, namely diagrams 3 and 4 of fig. ^ where one encoun- 
ters problems with the traces involving 75. There we meet the ambigous trace expression 

= S^'^^TrlS^^xlal^], (5.19) 

where S^^x is defined as before. All other ambigous expressions can be reduced to this one 
together with other non-ambigous terms. Besides the treatment of this problem in NDR 
as described in the previous section, we can attempt an Ansatz for Ua^ to replace it by 

U„ -96(l + t;£ + ...)7,75, (5.20) 

with V representing the ambiguity which is of 0{e). 

Although V does appear at intermediate steps of the calculation of diagrams 3 and 
4, after inclusion of evanescent contributions, it cancells out in the final expression. The 
result thus obtained agrees with the result obtained through the procedure of sect. 5.1. 
This gives us confidence that our treatment of NDR is correct. 

We would like to point out that the replacement of eq. (|5.20|) seems to be related to 
the procedure recently proposed by Korner et.al. We emphasize that in both works 



the scheme presented is a prescription, which although shown to work upto two-loops for 
certain processes, still needs to be proven consistent to all orders. 

5.3 Results 

(2) (2) ~ 
The singular terms in the diagrams Pi - with Q2 and Qq insertions in the NDR and 

HV scheme, are collected in tables ^|-^ The two-loop counter terms have been included 
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diagram by diagram. Colour factors have been omitted. The overall normalization is such 
that after the inclusion of colour factors separate contributions to (&2)ij of eq. ( p. 14 ) are 
obtained. The singularites in tables include also contributions from penguin diagrams 
obtained from fig. ^ by left-right reflections. It should be stressed that these additional 
diagrams give generally different singularities than the diagrams given explicitly in fig. ^. 
As usual, the result of a given penguin diagram contains two Dirac structures, 



V, 



LR 



7/. (1-75)^7'' (1 + 75) , 



V, 



LL 



7/.(l-75)®7''(l-75) 



(5.21) 



In diagrams P| 



(2) 



(2) 

PiQ the singular factors multiplying Vlr and Vll are the same, but 
(2) _ p>(2) observe, that diagrams Pi^\ Pg^^ 



, and PiQ do 



they differ in the case of Pii - Pi^ 
not contribute to the anomalous dimensions. 

The singular terms for Q2 and given in tables are also valid for Qi and Q5 
respectively. However, since Qi and Q5 are in the colour non-singlet forms, whereas Q2 
and Qq are in the singlet form, the colour factors differ. In fact, it turns out that in the 
case of the insertion of Qi and only the colour factors in diagrams 3, 4, 11, and 12 are 
non-vanishing, and only these diagrams contribute to '^^\Qi) and 7i^-*(Q5)- 

Including colour factors, using the results of tables and eq. (|2.14| ) without 02 (taken 
already into account in the current-current contributions), we obtain the expressions for 





1 

p 




1 

p 


7i^HQ5); 


, and 

p 





from which the full matrix 



can be found by means of formulae of sect. 5.1. 
In the NDR scheme we have 



7 



(1) 



7i,NDR(<5l> 



7i,NDR(Q2) 



7fNDR(Q5) 



3N 



_2_ 7 

3iV' ~ 3' 

32 86 176 
27 ^ 27iV2' ~W 



56 
27 



178 
27iV2' 



230 
27iV' 



r 

3. 

122 
'27 



94 86 130 \ 

— H 

27A^2' 27 27A^/ 



11 2 11 



16 250 70 

A^ H , — 

27 27A^ 27 



+ 



74 



110 



27A^2' 27 



A^ 



(5.22) 
(5.23) 
(5.24) 



^)/,(5.25) 
27A^/ ^ ' 



where the vectors are in the (Qa, Q4, Q5, Q^) space. 
In the HV scheme we have 



3Ar-l,l,-3Ar+l,l 



56 
27 



86 



27A^2' 27 



110,, 140 

A^ - 



27A^ 



128 



58 



27A^2' 27 



38 148 \ 
A^+ ttzttt) 



27A^. 



(5.26) 
(5.27) 
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3N + , -, 3A^ 

3N' 3' 



- - )/ 

3A^' 3/ ' 



Using next eqs. 
and 



7s,NDR 



(5.28) 

27 27A^2' 27^' ' 27N' 27 ' 27Ar2' 27^' 27A^y ' '^^'^^'^ 
STTP - (|5.6| ) in conjunction with (|5.9| ) - ( ^.121) , we find the full matrices 



128 



94 



20,, 202 



38 86 110 
~27^ 27iV2' '27 



N 



158 



7s,HV 



We observe that for large the contributions from penguin diagrams to 7^^^'' grow at 
most as A^. Since a'^ ~ 0{1/N'^) the full contribution of penguin diagrams to the third 
term in eq. ( |1 . 1|) vanishes in the Large- limit. Yet for = 3 the penguin diagrams play 
a considerable role in the numerical values for 

6 Full Two— Loop Anomalous Dimension Matrix 7g^^) 
6.1 Basic Result of this Paper 

Adding the current-current and penguin contributions to 7s^^^ found in sections 4 and 5 



respectively, we obtain the complete two-loop matrices 7s ndr 7s hv- 

We first note that 48 elements of these matrices vanish in both schemes. These are 
given by 



where 



3, 
7, 



/s.NDR 



.,10, J = 1,2} 
J = 9, 10} ; 



/s.HV 



= 

u 

1,...,6, j = 7,...,10} 

9,10, J = 7, 8} 



(6.1) 



(6.2) 



The remaining 52 elements of each matrix are non-vanishing. These entries in the 
anomalous dimension matrix are given in table | for arbitrary number of colours and 
flavours. 

For phenomenological applications we need only the results with A^ = 3. We give 
them in appendix C for an arbitrary number of flavours. It is evident from this table 
that the two-loop O^a"^) corrections to the anomalous dimensions of the operators Qi are 
substantial although they are quite different for the two schemes considered. 



6.2 Compatibility of NDR and HV Results 

Using the results of table 5 and the one- loop results for Afg of section 3 we can test 
whether the compatibility relation ( |2.15| ) is satisfled. This turns out to be indeed the case 
which constitutes a test of our calculation. 
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6.3 An Additional Test of the Calculation 



On general grounds, renormalizability relates the singularities of degree k in the n-th. order 
of perturbation theory with the singularities of degree [k + 1) in the (n+ l)-th order. This 
fact allows for an additional check of our computation by considering the corresponding 
relation between the coefficients of the singularities at two-loops and the 1/e singular- 
ities at one-loop. Demanding the finiteness of the renormalized Green function T^^\ and 
the anomalous dimension matrix 7, a procedure analogous to the derivation of eqs. ( p.l3| ) 
and ( |2.14| ) leads to 

f 22 = - ^ 7i°^' +{\f3o- ^1)7^°^ + [i/3o + 5o) a,-2al]i, (6.3) 

where f22 is the coefficient of the 1/5^ singularity at O^a"^) of the bare Greens function 
and 60 is the leading order coefficient of the renormalization group function associated 
with the gauge parameter. In the Feynman gauge it is given by 

60 = -^N+p. (6.4) 

The values for Pq and ai were already given in eqs. ( |2.18D and ( |2.19| ) respectively. 

However, to perform the test we have to add the contributions of the diagrams of fig.^ 
to the result which can be obtained by means of tables Despite the fact that they 
have no 1/e divergence and hence do not contribute to the anomalous dimension matrix, 
they contribute to f 22. Their contribution to f 22 is given by 

Af22 = ^(/?o + Cp)7f , (6.5) 

where the term proportional to (3o stems from the renormalization of the gluon-propagator 
and the term proportional to Cp from the vertex renormalization. Let us recall that we 
work in the background field gauge, where the /3-function can be obtained directly from 
the renormalization of the gluon propagator ||30[| . 

Putting everything together, we find that the relation of eq. (|6.3| ) is indeed satisfied. 
In addition, similar to the test of the compatibility of different schemes, this relation could 
also be split in a part solely originating from current-current diagrams which is satisfied 
separately, and a relation involving mixing between current-current and penguin diagrams. 
We will however not elaborate on this any further. 



7 Summary 

We have presented the details and the explicit results of the calculation of the two- 
loop anomalous dimension matrix 0{ag) involving current-current, QCD-penguin and 
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electroweak-penguin operators. Performing the calculation in two schemes for 75, NDR 
and HV, we have verified the compatibility of the anomalous dimension matrices obtained 
in these two schemes. We have shown how the use of two different operator bases allows 
to avoid a direct calculation of penguin diagrams with closed fermion loops. This enabled 
us to find in an unambiguous way the two-loop anomalous dimensions in the simplest 
dimensional regularization scheme, the one with anticommuting 75. The results of this 
paper generalize our earlier paper |]14| where electroweak penguin operators have not been 
taken into account. The two-loop anomalous dimensions matrix presented for the first 
time in ref. fl^ and generalized in section 6 of the present paper will play a central role in 
any analysis of non-leptonic decays of hadrons which goes beyond the leading logarithmic 
approximation. A consistent next-to-leading order analysis involving electroweak penguin 
operators requires the calculation of ^^g-*, i.e. the two-loop anomalous dimension matrix 
0{aas). A detailed account of this calculation is given in |jl5[. The full renormalization 
group analysis with the anomalous dimension matrix of eq. ( p..l|) will be presented soon in 
ref. [|l6l. There the numerical values for the Wilson coefficient functions Cj of the operators 
Qi including next-to-leading order corrections can be found. 
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Appendices 

A One— Loop Anomalous Dimension Matrix 7g^'^) 
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B Tables of Singularities and Two— Loop Anomalous 
Dimension Matrix 7g^^) 

Tables list the 1 / e"^- and 1 / ^-singularities of the various Feynman diagrams obtained 
in the NDR and HV scheme. The values given already include diagram by diagram two- 
loop counterterms. Colour factors are omitted and a common overall factor a'^/{47r)'^ is to 
be understood. 

In table |I| the column labeled D refers to the numbering of current-current diagrams 
in figs. ^ and ^. Similarly, the P-column in tables ^|-^ refers to penguin diagrams of figs. ^. 
In addition, in table ^ the column labeled M gives the multiplicity of each diagram, which 
already has been included in the singularities quoted. The singularites in tables include 
also contributions from diagrams obtained from fig. ^ by left-right reflections. 
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Table 5: Full QCD Anomalous Dimension Matrix for the NDR and HV scheme 

(with vanishing entries omitted tacitly). 
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Table 5: Pull QCD Anomalous Dimension Matrix {'yP)ij for the NDR and HV scheme 
(continued; with vanishing entries omitted tacitly). 
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C Two— Loop QCD Anomalous Dimension Matrix 7g^^^ 
in NDR and HV Schemes for TV = 3 
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D Figures of Feynman Diagrams 




(1) (2) 

Figure 2: 
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Figure 5: 
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Figure 6: 




Figure 7: 
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Figure Captions 



Figure 1: The three basic ways of inserting a given operator into a four-point function: 
(a) current-current-, (b) type 1 penguin-, (c) type 2 penguin-insertion. The curled 
hues denote gluons. The 4-vertices "0 ®" denote standard operator insertions. 

Figure 2: One-loop current-current diagrams contributing to 7^^°^. The meaning of lines 
and vertices is the same as in fig. 1. Possible left-light or up-down reflected diagrams 
are not shown. 

Figure 3: One-loop type 1 and 2 penguin diagrams contributing to The meaning of 
lines and vertices is the same as in fig. 1. 

Figure 4: Two-loop current-current diagrams contributing to The meaning of lines 
and vertices is the same as in fig. 1. In addition shaded blobs stand for self-energy 
insertions. Possible left-light or up-down reflected diagrams are not shown. 

Figure 5: Two-loop penguin diagrams contributing to The curled lines denote glu- 
ons. Square-vertices stand for type 1 and 2 penguin insertions as of figs. |l](b) and 
(c), respectively. Possible left- light reflected diagrams are not shown. 

Figure 6: Two-loop penguin diagrams having no 1/e divergence and hence do not con- 
tribute to the anomalous dimension matrix 7^^^^ The meaning of lines and vertices 
is the same as in fig. 5. In addition shaded blobs stand for self-energy insertions or 
vertex corrections. 

Figure 7: Two-loop penguin diagrams vanishing identically in dimensional regular izat ion. 
The meaning of lines and vertices is the same as in fig. 5. In addition shaded blobs 
stand for self-energy insertions or vertex corrections. 
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